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cn ' Abstract. A brief outline is given of the description of phase transition 

^ ■ kinetics in condensed matter systems with a continuous symmetry, empha- 

t:;j- . sising the roles of dissipation, coarse-graining and scaling. The possible 

^ ! relevance of these ideas to the early universe is explored in the contexts of 

1-^ ' the GUT string transition and the electroweak transition. 

^, ■ 

•^ : 1. INTRODUCTION 

c^ I How fast do phase transitions occur? Remarkably, it is found that although it is often 

straightforward to estimate a characteristic relaxation time for the microscopic degrees 
of freedom, the actual characteristic time for completion of the phase transition may 
be many orders of magnitude greater. For example, laboratory experiments indicate 
that following a temperature quench, the transition to the superconducting state of 
a normal metal in a magnetic field may take many minutes. The primary reason for 
the slowness of the transition is the formation, interaction and subsequent dynamics of 
topological defects. In systems with a discrete symmetry and a scalar order parameter, 
such as binary alloys, the topological defects are domain walls, whilst in systems with 
a continuous symmetry and a vector or tensor order parameter, such as certain liquid 
crystals, the defects may be strings and monopoles. The motion and mutual annihilation 
of the defects is usually the rate-determining step for the transition, and is affected by 
such factors as the presence of dissipation or disorder, the range of the interactions 
between defects, and even their homotopy classification. 
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Such phenomena are ubiquitous in condensed matter, yet only recently have de- 
tailed studies been made of phase transition kinetics in condensed matter systems with 
non-scalar order parameters. Although interesting in their own right, these exam- 
ples may be regarded as caricatures of the phase transitions that are thought to have 
occurred in the early universe as it cooled. The principal analogous feature is the spon- 
taneously broken global continuous symmetry and, in superconductors at least, the 
existence of a local gauge symmetry. 

The purpose of this article is to outline briefly the way in which the kinetics may 
be investigated theoretically, with an eye towards using the techniques of condensed 
matter physics in a cosmological context, where appropriate. In particular, the work 
of the group at Illinois on phase ordering in systems with continuous symmetries is 
relevant. This primarily numerical work includes studies of the non-conserved dynamics 
of the XY model in two dimensions fll], three dimensions [§], the conserved XY model 
in one and two dimensions 0, the dynamics of the superconducting transition |Q, the 
dynamics of the Ising gauge theory, where there is no local order parameter and 
ordering in uni- and bi-axial liquid crystals^. Other relevant numerical studies are 
those of Toyoki0 and Bray, who has recently given a complete review of the topic of 
phase ordering in systems with continuous symmetries 0. 

Since all of the literature and an extensive review are easily available, it does 
not seem worthwhile to provide duplication of the results here. Instead, the focus will 
be on principles and concepts, together with some remarks comparing the procedures 
used in condensed matter physics with those used in cosmological applications. These 
are particularly pertinent in the case of the electroweak transition, where the phase 
transition kinetics may not occur via nucleation and growth, as has sometimes been 
assumed. 

1.1. Scaling 

Let us consider qualitatively the sequence of different time regimes exhibited by a binary 
alloy undergoing spinodal decomposition. First, the small amplitude long- wavelength 
fluctuations of the order parameter, present from the initial conditions, become ampli- 
fied exponentially: this behaviour is predicted by the linearised equations of motion, 
and is rarely observed. The nonlinearities in the equation of motion quickly stabilise 
the order parameter at its two equilibrium values almost everywhere. Thus a series 
of domains has been formed, separated by domain walls. This interlocking pattern of 
domains subsequently coarsens, driven by the excess energy from the curvature of the 
domain walls. At very long times, the system attains equilibrium, which energetically 
should simply be a single domain wall dividing the system into two coexisting equilib- 
rium phases. The intermediate time regime, where the domains simply coarsen, exhibits 
dynamic scaling^. 

Dynamic scaling simply means that at large enough times t, the emerging pattern 
contains only one time-dependent length scale L{t). Thus, the equal time two-point 
correlation function of the order parameter ip defined by 

C(r,t) = (^(r,t)^(0,t)) (1) 

is actually only a function of a reduced variable: 

C(r,t)=F(r/L(t)) (2) 



where F is known as a scaling function. In practice, for finite times, there is a weak 
exphcit time dependence in F, which proceeds to increasingly short distances, whilst 
L{t) converges to its 'ideal' power law form, as discussed below. 

Of course, there are other length scales present: the bulk correlation length for 
longitudinal order parameter fluctuations ^(T), which is a function of the temperature 
T after the quench, and the microscopic length scale a, on the order of molecular dimen- 
sions. But for sufficiently long times, the following inequalities hold: L{t) ^ ^(T) > a. 
Since thermal fluctuations are operative up to the scale of ^, the unimportance of ^ 
relative to L is sometimes stated in a suggestive way by saying that phase ordering "is 
controlled by a zero temperature fixed point" . (Despite the intuition that some sort of 
renormalisation group (RG) approach should form the basis of a general theory for the 
approach to equilibrium, a predictive RG theory remains elusive.) The scaling regime 
may be called an intermediate asymptotic regime. For example, in a binary alloy, at 
short times, the initial fiuctuations in the system are being amplified by the unstable 
growth process, and domain walls are formed, whereas for long enough times in a finite 
system, the system reaches thermal equilibrium, where only one domain wall traverses 
the system. In both of these regimes, thermal fiuctuations are important. For some 
set of intermediate times, whose duration is an increasing function of the system size, 
the scaling regime is observed, and thermal fiuctuations are not important in the sense 
described above. Empirically, it is found that the growth of the characteristic scale L 
follows a power law form: 

m = A^T) {t/r{T))\ (3) 

where t{T) is the temperature dependent order parameter relaxation time, A is an 
amplitude assumed to be of order unity (but not yet measured as far as I know), 
believed to be universal and cj) is now believed to have the value 1/3 for alloys or other 
conserved systems with a discrete symmetry in dimension d > 2. For systems with a 
non-conserved parameter in these dimensions, has the value 1/2. Eq. |^ is valid near 
the transition temperature Tc, and a more general form valid for all T < T^ has been 



proposed by Bray [10 



Heuristic arguments, given elsewhere in this volume (see, for example, the lectures 
by Bray), connect these power laws with the motion of domain walls. So one is led to ask 
how these scaling results are affected when the system exhibits a continuous symmetry, 
rather than a discrete one, so that the domain walls no longer exist. In fact, other 
topological defects are present in such systems, and can give rise to different growth 
laws. These considerations formed the motivation for the work done at Illinois. 

1.2. Computer simulation 

In order to observe quantifiably the scaling regime in a computer simulation, two con- 
ditions must be met: the system size needs to be as large as possible, so that the regime 
lasts as long as possible, and then the longest possible times must be attained. In addi- 
tion to these requirements, a number of other results are known which are diagnostics of 
the scaling regime: that is, they are only satisfied in the scaling regime. These are the 



Tomita sum rule|Tl|], Yeung's law for the A; ^ behaviour of the X-ray scattering form 



factor at small wavenumbers A;|T^, and Porod's law for the form factor's short distance 
behaviour ||13||. In the absence of an obvious small parameter, systematic analytical 
work has rarely been possible — a notable exception is the study of the 0{N) model 
for iV — >■ CX3, where there are no topological defects fl^ — and much of our knowledge 



has come from computer simulations. Thus, in order to estimate accurately (e.g.) the 
exponent 0, computationally efficient techniques have been designed to probe as far into 
the scaling regime as possible: these are discussed in section 3.3. In fact, rather little 
of the numerical work to date actually satisfies the known criteria for the asymptotic 
regime; the most complete work to date on the alloy phase separation problem in three 
dimensions is that of Shinozaki and OonoHTSl. 



2. FORMULATION 

2.1. Level of description 

The main difficulty in constructing a theory for phase ordering is the disparity between 
the different length scales present (^, L{t) and a), and the complexity of the actual 
microscopic equations of motion. However, we are mainly interested in phenomena on 
the scale of L, so a coarse-grained description is adequate. This sort of approach is 
common in condensed matter physics: for example, the BCS theory of superconductiv- 
ity is a well-tested microscopic theory but is virtually useless in situations with spatial 
variation, such as near boundaries or for time dependent phenomena. Instead, the 
Ginzburg-Landau theory, a phenomenological theory for coarse-grained order parame- 
ter, slowly varying on the scale of ^, is used. 

2.2. Coarse-grained order parameter 

Let us be more precise, taking as our example, the case of a superconductor. The order 
parameter \& is zero for T > T^ and nonzero for T < Tc. It can be defined in terms of 
an anomalous Green function by 



vl/(x)oc(^i(x)^T(x)), (4) 

where ipi{^) is a down-spin electron field operator and the angle brackets denote an 



equilibrium thermal expectation value |]T6|. The normalisation of \l/ may be chosen with 
a convention that does not concern us here. In the Meissner and normal phases of 
a superconductor, the system is translationally invariant and \1/ is spatially uniform. 
Near equilibrium, ^ may vary in space, but we will only consider the long wavelength 
variations: in spinodal decomposition, the instabilities occur at wavelengths long com- 
pared with a. Conceptually, the coarse-grained order parameter is defined by the long 
wavelength Fourier components '^\f_: 

v1/a(x) = Y. e^'^'^k (5) 

|k|<A 

where the coarse-graining scale A satisfies a -C A^^ < ^. We want to coarse- grain 
as much as possible, but if we coarse-grain beyond the correlation length then we will 
have two-phase coexistence within one coarse-graining volume. This over-coarse-grained 
description would then not be able to tell us about the dynamics of phase separation. 
The coarse-grained free energy governs the effective dynamics of these long wavelength 
modes, and is obtained by integrating out the short wavelength modes up to the scale 
A: 

g-FA{*A(x)}ABT ^ f Yl (i^ke-^/'^^(5(^A(x) - Y. e^^'^k), (6) 

|k|>A |k|<A 



where H is the Hamihonian and ks is Boltzmann's constant. Thus, for a given coarse- 
grained order parameter profile in space, a value for the coarse-grained free energy can 
be calculated. In thermal equilibrium, the probability that the coarse-grained order 
parameter has a given profile \1/a(x) is proportional to exp (— FA{\l/A(x)}/fcBT). The 
explicit calculation of Fa is technically complicated, and unnecessary to perform if one 
is interested in universal equilibrium quantities. If one is interested in computing the 
dynamics, certain aspects may not be universal, and will depend upon Fa. A recent 
example where Fa has been explicitly calculated using the RG is the work of Alford 
and March- Russell |1^ . 



It is important to notice the distinction between the calculation of Fa and the 
finite temperature effective action S of field theory. The latter involves no notion of 
coarse-graining {i.e. A = oo), and is rigorously convex. On the other hand, there is 
no requirement of convexity for the coarse-grained free energy, and indeed, it has the 
familiar double- well or wine-bottle form for T <Tc. The non-convexity of the one- loop 
effective action is an artifact of the loop expansion. 

2.3. Dynamics 

How should we describe the dynamics of the phase transition? Ideally one would solve 
for the time- dependent density matrix, using the exact microscopic Hamiltonian. Ob- 
servables such as the coarse-grained order parameter, or its correlation functions, could 
then be computed. This would involve averaging over the initial conditions with the 
appropriate Boltzmann weight, and modelling the quench itself. In other words, one 
would compute the exact dynamics, then coarse-grain. An alternative procedure would 
be to find the effective equations of motion for the coarse-grained order parameter by 
directly coarse-graining the equations of motion for the density matrix. A plausible 
substitute is to write down the phenomenological equation of motion for \&a5 assuming 
that the driving force is Fa. Even this is rarely done, because usually Fa is not known. 
Instead, a phenomenological form of Fa is used. 

A phenomenological Langevin equation for \1/a is obtained by assuming that the 
force driving the system towards thermal equilibrium is proportional to the deviation 
from equilibrium: 

^»^'*^ = -jt + '' P> 

where we have assumed that \1/a is a complex scalar field (appropriate for a superfluid, 
for example), and included a thermal noise term r] so that the system is guaranteed 
eventually to attain the global minimum of Fa. We shall refer to eq. |^ as the time- 
dependent Ginzburg-Landau equation (TDGL). As mentioned above, the noise term is 
believed not to play a significant role during the scaling regime. The relaxation time 
Tq is in principle calculable from a microscopic theory (such as BCS theory in the case 
of a superconductor). The Langevin equation above does not include any conservation 
laws, although these are easy to include when required. For a superconductor, an appro- 
priately generalised local gauge invariant time-dependent Ginzburg-Landau description 
can be written down (see e.g., [Q). 

It is difficult to assess the regime of validity of eq. 0. Usually it is regarded as 
a minimal model of phase ordering kinetics, in the sense that it correctly predicts all 
universal phenomena. However, it may not be accurate for other quantities, and in this 
sense, it is sometimes said that eq. |^ is a semiquantitative description. An important 



point is that we have assumed that the dynamics is purely relaxationaL Whilst this 
is expected to be valid near a critical point, where the correlation length is large, for 
quenches to low temperatures this is not necessarily the case. 

Let us use again our example of the coarse-grained dynamics of a superconductor 
to see what can go wrong. There, the microscopic BCS theory may be formulated 
as a set of coupled equations for the Green's functions of the theory, which in turn 
can be reexpressed as an integral equation for the order parameter ||16||. The TDGL is 
obtained by expanding the integral equation in powers of wavenumber k and frequency 
u. However, only at low temperatures, near the critical temperature or when there 
is strong scattering from impurities does this procedure yield a result independent of 
the ratio u/k, leading to a local partial differential equation. In all other situations, 
there is no possible local description of even the long wavelength, low frequency order 
parameter dynamics. Furthermore, at low temperatures, it turns out that the dynamics 



is not overdamped but exhibits wave-like solutions [|T8|. 



3. NUMERICAL METHODS 

3.1. Monte Carlo Simulation 

This is probably the most straightforward conceptually. Here, the Monte Carlo time is 
assumed to be proportional to real time, and the dynamics that the system undergoes 
is assumed to be somehow similar to that of the real system. Of course, this may not be 
the case: Monte Carlo simulation is only a stochastic process which samples the correct 
equilibrium distribution. However, a judiciously chosen Monte Carlo dynamics may 
well be a reasonable caricature of the actual dynamics. This method is usually very 
slow, and it has not proved to be reliable in extracting the correct long time behaviour. 
It is not the method of choice nowadays. 

3.2. Molecular Dynamics for the Defects 

As discussed at length in this volume, topological defects are created during the phase 
transition, and equilibrium may only occur after all such defects have annihilated. In 
many situations, the potential for interactions between defects may be calculated, and 
the defects treated as classical particles, subject to damping and the inter-defect po- 
tential. Accordingly the equation of motion for the particles may be readily solved. 



An example, in the case of superconductor dynamics, is ref. |19]. This approach is 
relatively fast in terms of computer time, and large systems may be straightforwardly 
treated. The disadvantage, however, is that it is usually non-trivial to obtain the cor- 
rect equation of motion for the defects, either because the potential may be hard to 
obtain, or because the damping may be difficult to calculate. An interesting example is 



the calculation of ordering in the non-conserved two dimensional XY model in ref. pO 
Here, the ordering proceeds via the annihilation of ± vortices, whose interaction poten- 
tial U is logarithmic in their separation R. The equation of motion for the separation 
R of two vortices is assumed to be 

dR dU 

"^ = -dR (') 

where z/ is a damping coefficient. This calculation, as well as the molecular dynamics 
simulation with many vortices, finds that the average separation L varies as L ~ t^/^. 



whereas the correct result is now beheved to be L ~ (t/logt)^/^: the logarithmic factor 
arises because the damping coefficient v is actually logarithmically dependent on R 



21| . Note that the above estimate of scaling exponents, based upon the overdamped 
equation of motion for a pair of topological defects, works well in other situations, even 
correctly predicting the crossover in a superconductor with penetration depth A from 
L ~ (i/logt)i/2 (L < A) to L ~ logt (L > A)|]. 

3.3. Order Parameter Evolution - PDEs and CDS 

Probably the most effective way to explore phase transition kinetics is to solve directly 
the TDGL or equivalent partial differential equation (PDE). This involves discretising 
the PDE on a space-time lattice: 

^(x, t) -^ ^(nAx, iAt) = ¥^ (9) 

and using a time-stepping algorithm such as the explicit Euler scheme 

6F 



¥+^ = ^1 + At 



5^* 



(10) 



Of course, eventually the continuum limit must be taken: Ax, At — > 0, subject to 
possible stability criteria giving an upper bound to At/Ax^. The disadvantage of this 
method is primarily that of speed and memory: to explore asymptotically large times, 
with a time step tending to zero requires many iterations of eq. |T0|. 

The cell dynamic system (CDS) method [p^ exploits universality in order to over- 
come this problem, and is now the most widely used approach in studying phase or- 
dering and other pattern formation problems in condensed matter. The basic idea is 
the observation that it is rather wasteful to model the phenomenon in question by a 
PDE, which one must then discretise to obtain a set of coupled maps (such as eq. |10D 
suitable for numerical computation. Instead, the phenomenon is modelled directly by 
a set of coupled maps, defined on a coarse-grained space-time lattice, with spatial cells 
of dimension of order 27rA~^. 

To illustrate the basic idea, consider the case of an order- disorder transition, where 
the order parameter is a simple non-conserved real scalar field obeying the TDGL, and 
equilibrating with the coarse-grained free energy 



F{*(x)}= fd'^x 



2^ ^2 4 



(11) 



The parameter a changes sign at T^, is initially positive when T > T^, and at t = 0+ 
is supposed to become instantaneously negative. When substituted into the TDGL, 
a nonlinear PDE is obtained. The CDS approach need make reference neither to the 



TDGL nor to the coarse-grained free energy functional of eq. |TT|. Let us consider the 
dynamics of the order parameter in one cell, throughout which the order parameter 
value is essentially constant, but allowed to vary in time. The time dependence of \E' 
will have three fixed points, two sy mme trically placed about one at 0, corresponding 
to the three extrema at \E' = 0, ±Ja/b. The former is an unstable fixed point of the 
cell dynamics, whereas the other two are stable fixed points. We can model the cell 
dynamics by any map A^At that has this fixed point structure. Thus 

^L^' = MAtiK}- (12) 



It is often convenient to take the map 

A^At{^} = l-3tanh(^), (13) 

although even simpler piece-wise linear maps have also been used. If desired, one could 
explicitly calculate A^a* from the TDGL, although there would be no point for present 
purposes: the phenomenological differential equation or phenomenological PDE have 
no more privileged status than our phenomenological map. The interaction between 
cells should reflect the role of diffusion processes, and should be as isotropic as possible. 
Since diffusion is simply a local averaging procedure, we couple the cells in the final 
form of our CDS by writing 



K^' = MAtiK] + D ((^y - ^'nj (14) 

where D is a phenomenological coefiicent, akin to a diffusion constant (but incorporating 
the spatial and time discretisation units) and the averaging operation is in {e.g.) two 
dimensions 

^ n.n. ^^ n.n.n. 

Here n.n. means nearest neighbours, n.n.n. means next nearest neighbours. This form 



of Laplacian is more isotropic than the conventional discretisation |]T5|, ^, as can be 
seen by examining isocontours of its Fourier-transform. 

Empirically, this CDS modelling and its extensions works very well in a wide 
variety of situations, and gives results that in some cases have been compared with 
direct integration of the TDGL. This is not surprising: the actual form of the map used 
is of little consequence for large scale structure, and only influences the detailed form of 
the order parameter profile near a domain wall or topological defect. Of course, since 
the free energy of eq. |Tl] is only phenomenological, its predictions for this variation are 
no more reliable than those of the CDS map. 

The key point about the CDS method, and one that is frequently misunderstood, 
is that there is no continuum limit: Ax and At are not infinitesimals. Thus, very rapid 
simulations are possible on large systems. The correct way to view the CDS map is 
that one has integrated or coarse-grained some microscopic equations of motion up to 
the space-time scale of interest, thus obtaining the map we have guessed phenomeno- 
logically. This is conceptually different from the usual approach of numerical analysis, 
which is to sample a PDE on a sequence of finer and finer meshes. There is no proof 
that the CDS algorithm is in the same universality class of the PDE (whatever that 
may be), but the important point is that both are in the same universality class as 
the physical phenomenon of interest. This is the job of the physicist: to identify and 



characterise such universality classes. Ref. p4| describes recent work on the application 
of RG to extract universal features of PDEs. 

In the literature, a popular variant on the CDS method is often encountered, in 
which the conventional Euler discretisation of a PDE is used, but with large At and Ax. 
This works well too, although one must then be careful not to identify the potential 
used in this discretisation of the PDE with the potential in the original continuum limit 
PDE, for the reasons described in the preceding paragraph. The tanh map advocated 
above is simply a convenient form of the potential, which happens to avoid a secondary 
numerical instability by virtue of being injective. 



4. COSMOLOGICAL PHASE TRANSITIONS 

In this section, I will present some remarks and observations on several issues of interest 
to cosmologists. These arose during discussions at the workshop, whilst presenting the 
procedures and insights obtained from simulating string formation, evolution and the 
dynamics of local gauge theories in condensed matter, and from inspection of computer 
animations of phase transition kinetics with non-scalar order parameters. 

4.1. Violations of scaling 

Although scaling is observed generically in the presence of a conservation law, the 
nonconserved order parameter dynamics only leads to scaling solutions for symmetric 
or critical quenches. Thus, if the order parameter obeys the TDGL with the coarse- 
grained free energy of eq. |Tl|, then a symmetric or critical quench is one where the 
spatial average (^l/(x, 0)) = initially. However, it is not obvious that in a cosmological 
context, the ordering field will always satisfy this condition. For example, in a GUT 
transition at which strings are formed, earlier symmetry breaking transitions may create 
a bias in the Higgs field. In the toy-model case of the three dimensional XY model, 
simulations[0 confirm the prediction from mean field theory [p5| that scaling is violated, 
with the initial string distribution rapidly breaking up into small loops. The total length 
of string £[t) a time t after the quench is given by 

£(t)~riexp(-62t3/2) (16) 

where the bias is given by 6^ = | (^(x, 0)) p. For times less than the crossover time 
tc ~ b~^^^, the dynamics of the string network resembles that of the critical quench, but 
for longer times, it becomes readily apparent that there are two length scales present: 
the constant mean size of small loops, and the increasing mean distance between them. 
The success of the analytic mean field theory is actually typical. Although analytic 
approximations on the TDGL itself have not been successful so far, transforming the 
equation into an effective equation for the defects has proven to be the appropriate 
starting point for semi- quantitatively accurate predictionspG . 



4.2. The Kibble mechanism in the GUT transition 

It comes as a shock to most condensed matter physicists to learn that cosmologists 
have a name for the mechanism in which topological defects are formed during the 
quench through a phase transition. (To my knowledge, it has only occurred to one 
condensed matter physicist that topological defects would not be created during a 
rapid quench!) During a putative phase transition in a grand unified theory (GUT), 
string-like topological defects can be formed, which, many expansion times later, have 
been proposed to act as seeds for galaxy formation. Although the initial string density 
does not turn out to be important for existing models of structure formation, because 
the string network is believed to attain a scaling solution at long times, some attention 
has been given to the question of determining this quantity. 

By the term "initial string density" is meant the length of string per unit volume 
"just after the quench" . This only has a unique meaning if the quench can be considered 
to be instantaneous: otherwise, the order parameter configuration is determined by the 
history of the system, and there is no unique value for this quantity. 



Computer simulations and experiments on liquid crystals very clearly show that 
following an instantaneous quench from above T^, the order parameter fluctuates wildly 
with position, and the order parameter configuration is so disordered that it is meaning- 
less to identify topological defects. After all, a topological defect is a defect configuration 
in a smooth, ordered background, and when this does not exist, individual defects can- 
not be distinguished. This is equivalent to saying that the cores of the defects overlap. 
Thus, right after the quench, the initial string density pi is not well-defined. In fact, 
it first becomes well- definable at a time ti when the average separation R{ti) between 
defects is greater than the core size, given by the order parameter correlation length 
^{Tf) at the final temperature Tf reached by the quench. Hence, the "initial" string 
density (meaning the string density at this time) can only be given by 

p. = l/R{Ufr^i{Tf)-\ (17) 

Note that the string correlation function at this and subsequent times, which will depend 
on the quench history, as well as the static correlations in the order parameter from the 
starting temperature above Tc- 

4.3. Electroweak transition 

There is currently considerable interest in the dynamics of the electroweak transition, 
due to the suggestion that baryon asymmetry arises in the propagating bubble walls 
accompanying this putative first order transition. For a review, see the article by Turok 
in this volume. The electroweak transition is conceptually similar to the superconduct- 
ing transition when gauge field fluctuations are included [^, and we shall couch our 
discussion in these terms. 

The equilibrium statistical mechanics near the superconducting transition is ob- 
tained from the partition function 

where the coarse-grained free energy is usually taken to be 

F{^,A}= /d'^x [V^I' + al^l' + Vl' + ^^^^^l- (19) 

J 2 2 Stt 

Here A is the electromagnetic vector potential, D"^ = V — ieA is the covariant deriva- 
tive, and we have worked in the gauge where the scalar potential is zero and V ■ A = 0. 



Halperin, Lubensky and Map^ showed that the functional integral over the gauge field 



leads to a term in the resultant coarse-grained free energy density 

F{^} = -keT log f DA exp(-F{$, A}/A;bT) (20) 

which is proportional to — I^E'P, and thus generates a first order transition, which gener- 
ically proceeds by nucleation and growth of bubbles. For strong supercooling, the 
nucleation barrier vanishes, and the dynamics proceeds by spinodal decomposition. In 
the cosmological context of the electroweak transition, it is generally believed that the 
transition occurs via nucleation and growth. This occurs because the Higgs mechanism 
causes the Higgs field to lose all degrees of freedom apart from the modulus |\1/|, and 
it is this quantity which appears in the effective potential, and which varies across a 



domain wall. There are well-known problems associated with the computation of the 
effective potential in this context, and the nature of the transition is currently unclear. 
The main remark that I would like to make here concerns the dynamics of the 
phase transition. In general, it is not correct to use the derivative of F{'$} as the 
driving force in a TDGL for "$. The reason is that the dynamics of \E' depends upon 
the dynamics of the gauge field, and at least in the case of superconductors, the rel- 
evant time scales for relaxation of both \E' and A are of the same order. The correct 
procedure is to coarse-grain both the microscopic order parameter and gauge field, and 
then to solve the coupled dynamical equations for both \E' and A. In the case of the 
superconductor transition, this has been done in ref. [ffl with the result that growth 



does not proceed by nucleation and growth of bubbles. Elder and myself |]2^ have re- 
peated the calculation using the effective potential F, in which the gauge field has been 
completely integrated out, and checked that the results are quite different from those 
obtained using the full dynamical equations for both the order parameter and gauge 



field. A similar result was obtained by Ye and Brandenberger|^9| in their numerical 
simulations of the Abelian Higgs model, where no evidence for domain wall formation 
was observed. If this conclusion is supported by further investigation, then there could 
be important ramifications for proposed mechanisms of electroweak baryogenesis. 

In conclusion, it seems that the time is ripe for a systematic study of phase 
transition kinetics and the dynamics of fields, not only in the context of condensed 
matter, but also in other areas where the space-time behaviour of fields is of interest. 
These include cosmology and perhaps the physics of the quark-gluon plasma, soon to 
be probed by RHIC. 
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